FEXPMaterial.h

Go to the documentation of this file.

// © FEXP, FEXPEnterprise Solver, Ing. Vaclav Rek
// Library for solid materials.
// Compiler must support C++ ver.14 and later
#ifndef _CFEXPMATERIAL_H_
#define _CFEXPMATERIAL_H_
#include "FEXPCommon.h"
#include "FEXPDataContainer.h"


enum EMaterial
{
  eNoMaterial = 0,
  eSimpleElastic
};

#define WAVE_SPEED_LAMBDA(statm_lambda) [](t_fexpcommon_ct E, t_fexpcommon_ct rho, t_fexpcommon_ct nu) statm_lambda
#define WAVE_SPEED_PAIRMP(key, statm_lambda) { key, WAVE_SPEED_LAMBDA(statm_lambda) }

class ICFEXPMaterialDynCharacter
{
public:
  virtual ~ICFEXPMaterialDynCharacter() { }

  // pure virtual member functions
  virtual t_fexpcommon_ct GetWaveSpeed(FEXPFEDefinition::EFEXPFiniteElementGeom egtype) = 0;
  virtual t_fexpcommon_ct GetMatDnsRho()                                                = 0;
  virtual t_fexpcommon_ct GetMatYngMod()                                                = 0;
protected:
  // wave base speed dynamic characteristics for respective FE type
  std::map<FEXPFEDefinition::EFEXPFiniteElementGeom, std::function<t_fexpcommon_ct(t_fexpcommon_ct, t_fexpcommon_ct, t_fexpcommon_ct)>> _wave_base_characteristics =
  {
    WAVE_SPEED_PAIRMP(FEXPFEDefinition::EFEXPFiniteElementGeom::eFE_1D_in2D, { return sqrt(E / rho);                                                  }),
    WAVE_SPEED_PAIRMP(FEXPFEDefinition::EFEXPFiniteElementGeom::eFE_2D_in2D, { return sqrt(E / (rho * (1.0 - pow(nu, 2))));                           }),
    WAVE_SPEED_PAIRMP(FEXPFEDefinition::EFEXPFiniteElementGeom::eFE_2D_in3D, { return sqrt(E / (rho * (1.0 - pow(nu, 2))));                           }),
    WAVE_SPEED_PAIRMP(FEXPFEDefinition::EFEXPFiniteElementGeom::eFE_3D     , { return sqrt((E * (1.0 - nu)) / (rho * (1.0 + nu) * (1.0 - 2.0 * nu))); })
  };
private:
  // [no private members  ] -----------------------------
};


class ICFEXPMaterialBase
  : public ICFEXPMaterialDynCharacter, public ICFEXPModelDataIntf
{
public:
  ICFEXPMaterialBase(size_t id) : _id(id) { }
  virtual ~ICFEXPMaterialBase() { }

  // pure virtual member functions for calculation
  virtual void CalculateStress(Ptr<CFEXPMatrix2D> strn, Ptr<CFEXPMatrix2D> d_strn, Ptr<CFEXPMatrix2D> strs, t_fexpcommon_ct dt) = 0;
  virtual Ptr<CFEXPMatrix2D> GetStress() = 0;
  // clone material for layers of layered FE
  virtual Ptr<ICFEXPMaterialBase> GetSameMaterial() = 0;
  // type of a solid material material
  virtual EMaterial GetMaterialId() = 0;

  // override pure virtual member functions
  virtual size_t GetId          ()          override { return _id;                                   }
  virtual size_t GetModelId     ()          override { return FEXPCOMMON_DEFAULT_VALUE;              }
  virtual size_t GetType        ()          override { return size_t(ESystemElementType::eMaterial); }
  virtual size_t GetThreadId    ()          override { return FEXPCOMMON_DEFAULT_VALUE;              }
  virtual size_t GetThreadNumber()          override { return FEXPCOMMON_DEFAULT_VALUE;              }
  virtual void   SetThreadId    (size_t id) override { }
protected:
  // [no protected members] -----------------------------
private:
  size_t _id;
};


class CFEXPSimpleElasticMaterial
  : public ICFEXPMaterialBase
{
public:
  CFEXPSimpleElasticMaterial(size_t id, t_fexpcommon_ct rho, t_fexpcommon_ct E, t_fexpcommon_ct nu);
  virtual ~CFEXPSimpleElasticMaterial() { }

  // dyn. char.
  virtual t_fexpcommon_ct GetWaveSpeed(FEXPFEDefinition::EFEXPFiniteElementGeom egtype) override;
  virtual t_fexpcommon_ct GetMatDnsRho()                                                override;
  virtual t_fexpcommon_ct GetMatYngMod()                                                override;
  // override compute
  virtual void CalculateStress(Ptr<CFEXPMatrix2D> strn, Ptr<CFEXPMatrix2D> d_strn, Ptr<CFEXPMatrix2D> strs, t_fexpcommon_ct dt) override;
  virtual Ptr<CFEXPMatrix2D> GetStress() override;
  // clone
  virtual Ptr<ICFEXPMaterialBase> GetSameMaterial() override;
  // type of a solid material material
  virtual EMaterial GetMaterialId() override { return ID; }

  // override pure virtual base functions for NETWORK serialization/deserialization
  virtual Ptr<t_SerializedData> NetSerialize  ()                           override { FEXPCOMMON_NOT_IMPLEMENTED_EX; }
  virtual void                  NetDeserialize(Ptr<t_SerializedData> data) override { FEXPCOMMON_NOT_IMPLEMENTED_EX; }

  static const EMaterial ID;
protected:
  // [no protected members] -----------------------------
private:
  void initialize();
  // material characteristics
  t_fexpcommon_ct _rho;
  t_fexpcommon_ct _E;
  t_fexpcommon_ct _nu;
  // material matrix
  static const size_t _dim_D;
  Ptr<CFEXPMatrix2D>  _material_D;
  Ptr<CFEXPMatrix2D>  _material_S;
};


class CFEXPRigidMaterial
  : public ICFEXPMaterialBase
{
public:
  CFEXPRigidMaterial(size_t id, t_fexpcommon_ct rho, t_fexpcommon_ct E, t_fexpcommon_ct nu);
  virtual ~CFEXPRigidMaterial() { }

  // dyn. char.
  virtual t_fexpcommon_ct GetWaveSpeed(FEXPFEDefinition::EFEXPFiniteElementGeom egtype) override;
  virtual t_fexpcommon_ct GetMatDnsRho()                                                override;
  virtual t_fexpcommon_ct GetMatYngMod()                                                override;
  // compute
  virtual void CalculateStress(Ptr<CFEXPMatrix2D> strn, Ptr<CFEXPMatrix2D> d_strn, Ptr<CFEXPMatrix2D> strs, t_fexpcommon_ct dt) override { }
  virtual Ptr<CFEXPMatrix2D> GetStress() override { return Ptr<CFEXPMatrix2D>(); }
  virtual Ptr<ICFEXPMaterialBase> GetSameMaterial() override;
  // type of a solid material
  virtual EMaterial GetMaterialId() override { return ID; }

  // override pure virtual base functions for NETWORK serialization/deserialization
  virtual Ptr<t_SerializedData> NetSerialize  ()                           override { FEXPCOMMON_NOT_IMPLEMENTED_EX; }
  virtual void                  NetDeserialize(Ptr<t_SerializedData> data) override { FEXPCOMMON_NOT_IMPLEMENTED_EX; }

  static const EMaterial ID;
protected:
  // [no protected members] -----------------------------
private:
  // material characteristics
  t_fexpcommon_ct _rho;
  t_fexpcommon_ct _E;
  t_fexpcommon_ct _nu;
};

#endif // !_CFEXPMATERIAL_H_